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Abstract. Particle production sources in pPb collisions at LHC energies are investigated in pseudorapidity 
space as functions of centrality. The nonequilibrium-statistical relativistic diffusion model (RDM) with 
three sources is applied to the analysis of preliminary charged-hadron distributions measured by the 
ATLAS and ALICE collaborations in pPb collisions at 5.02 TeV. The size of the midrapidity source 
relative to the fragmentation sources in pPb collisions is considered as function of centrality, and incident 
energy. At LHC energies particle production from the mid-rapidity gluon-gluon source is larger than 
that from the fragmentation sources at all centralities. Conclusions regarding the statistical equilibration 
process are drawn. 


PACS. 25.75.-q Relativistic heavy-ion collisions - 24.10. Jv Relativistic models - 24.60.-k Statistical theory 
and fluctuations 


1 Introduction 

Whereas charged-particle production in symmetric heavy- 
ion collisions has been investigated in great detail both in 
the AuAu system at RHIC c.m. energies of 19.6 GeV to 
200 GeV [T], and in PbPb at LHG energies of 2.76 TeV 
per particle pair [ 2 ], asymmetric systems such as pFh at 
the current LHC c.m.energy of 5.02 TeV are still calling 
for further exploration. 

So far, pseudorapidity distributions of produced charged 
hadrons in pFh have been measured in a rather limited 
range of pseudorapidities —2 < 77 < 1.9 (ALICE collab¬ 
oration 0) or 1 77 1 < 2.7 (ATLAS collaboration [3]), and 
centrality-dependent results are presently being assembled 

IHlS]- 

Here we investigate how these upcoming new results 
compare with the nonequilibrium-statistical relativistic dif¬ 
fusion model (RDM), which has already proven to be use¬ 
ful for detailed descriptions and predictions of charged- 
hadron distributions in symmetric systems at RHIC and 
LHC energies [ 6 ]. Its principal relation to other theoreti¬ 
cal models such as the statistical hadronization model 0 
IE], the color glass condensate ini, or hydrodynamics m 
has been discussed in a recent application of the RDM to 
symmetric systems [12] . The situation for the asymmetric 
case that is treated here is similar with respect to related 
models. 

In particular, the mid rapidity source in the RDM be¬ 
comes equal to the thermal (equilibrium-statistical) re¬ 
sult in the limit time to infinity. In the thermal model - 
which has a single source by construction - there is no 


limiting fragmentation at LHC energies m, whereas the 
RDM with two additional non equilibrium fragmentation 
sources does show limiting fragmentation [H] in agree¬ 
ment with the data. 

The longitudinal degrees of freedom are much further 
away from overall equilibrium than the transverse ones, 
which show equilibrium behaviour at small transverse mo¬ 
menta with collective expansion, and pQCD type events 
above pp — 7 GeV /c [15] . In this work, we integrate out 
the transverse motion, and consider only the longitudinal 
behaviour where the interplay of mid rapidity and frag¬ 
mentation sources for particle production becomes obvi¬ 
ous. 

In the RDM, the shapes of the distribution functions 
for asymmetric systems are much more sensitive to details 
of the model as compared to symmetric systems. Hence, 
valuable conclusions regarding the statistical equilibration 
processes in these relativistic systems can be drawn. More¬ 
over, the calculated results beyond the measured regions 
in (pseudo-)rapidity space provide detailed predictions for 
the forthcoming data analyses. 

The model is briefly explained in the next section, and 
in the main part of this note the results are compared with 
preliminary LHC data. 


2 The model 

The relativistic diffusion model (RDM) has been outlined 
in [16] , and applied to asymmetric systems at RHIC ener¬ 
gies in m, where dAu has been investigated at 200 GeV. 
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In the linearized version, the RDM is based on a Fokker- 
Planck type transport equation for the distribution func¬ 
tions Rk{y,t) in rapidity space, 


^Rk{y,t) = — \-^\{yeq-y)-Rkiy,t)\+D';-^Rk{y,t). 


dt 


dy L 




( 1 ) 



T1 


Fig. 1. (Color online) The predicted RDM pseudorapidity dis¬ 
tribution functions for charged hadrons in minimum bias pPb 
collisions at LHC c.m. energy of 5.02 TeV shown here are 
optimized in the mid-rapidity region to the preliminary AT¬ 
LAS (top) U], and published ALICE (bottom) data [3] (sys¬ 
tematic error bars only for ALICE). Corresponding y^-values 
per degrees of freedom are 0.436/45 for ATLAS, and 0.388/31 
for the ALICE data. The measured preliminary ATLAS data 
(0-90%) have been scaled by the geometric cross-section ratio 
cr(0 — 90%)/cr(0 — 100%) ~ 0.898 to make them directly com¬ 
parable to the minimum-bias ALICE results. The underlying 
distributions in the three-sources RDM are also shown. 


The total distribution function comprises an incoher¬ 
ent superposition of three sources: The source Ri repre¬ 
sents the fragmentation distribution arising from inter¬ 
actions of valence quarks in the projectile (proton) with 
\ow-x gluons in the target (Pb), R 2 the corresponding 
fragmentation distribution due to interactions of valence 
quarks in the target with gluons in the projectile, and 
Rs = Rgg the midrapidity distribution function arising 


from the mutual interaction of - mostly - low-x gluons in 
the fireball. The diffusion coefficients Dy account for the 
broadening of the distribution functions due to particle 
production and collisions, whereas the drift terms cause 
shifts of the distribution functions towards the equilib¬ 
rium value, and constrain the diffusion. The midrapid¬ 
ity source is already centered at the centrality-dependent 
equilibrium value of the rapidity, 


max ) + {m^{b)) ex.p{-ymax) 

2 ^ {ml{b)) ex.p{ymax) + {mj(b)) ex.p{-ymax) 

( 2 ) 

with the beam rapidities ^^eam = Tymax = T^^{^/SNN/rrip), 

the average transverse masses (^q2(^)) = 'y/^i,2(^) + 

and participant masses of the p- and Pb-like par¬ 

ticipants in pPb collisions that depend on the impact pa¬ 
rameter b. For very high energy such that exp {—ymax) — 0 
the equilibrium value of the rapidity becomes 


yeq{b) ^ 0.5 In 




( 3 ) 


For sufficiently large times t ^ 00 - in statistical equilib¬ 
rium - all three distribution functions tend to be centered 
at yeq{b). The relevant quantity that determines whether 
the centers of the fragmentation distributions reach y^q is 
the ratio of the centrality-dependent interaction time, and 
the source-dependent rapidity relaxation times, tinilTy^"^^ 
since for asymmetric systems 

(2/1,2W) = yeq[l-exp(-i/T^’^)]=F2/TOaa;exp(-i/T^’2). (4) 


It has turned out [6] that at RHIC and LHC energies, 
the interaction time in symmetric systems is too small for 
the system to reach statistical equilibrium [6] . It is only the 
central source that reaches equilibrium in rapidity space - 
although in transverse momentum space, the distribution 
functions of charged hadrons are actually very close to 
thermal equilibrium. The mean value of the moving mid¬ 
rapidity source remains {yl) = at all times, and 

the variance approaches equilibrium according to g\ (t) = 
DyTy[l — exp(—2t/T^)], for (5—function initial conditions 
as in [6]. 

The model has four parameters in case of symmetric 
systems: the ratio of interaction time and relaxation time 
(oi* the corresponding mean value as obtained 
from Eq. |4]), the variances FWHM A=2) of the 

fragmentation sources and (or FWHM Fgg ) of the mid 
rapidity source, plus the ratio of particles produced in the 
mid rapidity source, and in total for a given centrality. 
(If the total number of produced charged hadrons is not 
considered as experimental input, but as an additional pa¬ 
rameter, we have five parameters). These parameters are 
determined in minimizations of the analytical solu¬ 
tions with respect to the data. 

In case of asymmetric systems that are investigated 
in this work, the total distribution function is even more 
sensitive to the precise shape and the interplay of the 
three sub-distribution functions, together with the Jaco¬ 
bian transformation between rapidity and pseudorapidity 
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Table 1 . RDM parameters (mean values < yi ,2 > and widths G), and particle content of the sources for centrality-dependent 
pPb collisions at 5.02 TeV. See text for the number of RDM parameters. This y^-optimization (per numbers of degrees of freedom, 
dof) uses ROOT [18], preliminary ALICE charged-hadron dN/drj data [5], and equilibrium values yeq of the mid-rapidity source. 
The mean transverse momenta {pr) are taken from the data [19]. The beam rapidity is ybeam = Tl/max = T8-386. 


centrality 

(%) 

{pt) 

GeV/c 

Veq 

{yi) 

(j/2> 

A 

Ts 

^99 

iVi 

N2 

Ngg 

xVdof 

0-5 

0.800 

0.887 

-2.050 

2.900 

4.819 

4.500 

7.500 

100 

119 

308 

7.551 / 32 

5-10 

0.779 

0.886 

-2.090 

3.005 

4.591 

5.000 

7.994 

75 

104 

248 

7.967 / 33 

10-20 

0.767 

0.876 

-2.100 

3.100 

4.500 

5.100 

8.750 

53 

80 

236 

4.447 / 33 

20-40 

0.743 

0.818 

-2.300 

3.300 

4.353 

5.150 

8.900 

32 

47 

200 

1.461 / 32 

40-60 

0.713 

0.817 

-2.750 

4.500 

4.350 

6.000 

9.500 

18 

31 

157 

0.693 / 33 

60-80 

0.660 

0.563 

-2.760 

4.700 

4.200 

6.500 

10.000 

9 

18 

95 

0.463 / 33 

80-100 

0.608 

0.129 

-2.800 

4.710 

4.101 

6.998 

10.050 

4 

9 

31 

0.137 / 32 


space. As a particular example, the shoulders of the distri¬ 
bution functions are not symmetric: The fall-off depends 
on centrality, and is predicted to be steeper on the p-like 
side than on the Pb-like side, in particular in central colli¬ 
sions. Whereas this is indeed visible in dAu data at RHIC 
energies mini, LHC pPb data do not yet reach the re¬ 
gion in (pseudo-)rapidity space where this specific model 
prediction can be tested. 

For asymmetric systems, the mean value of the mid 
rapidity distribution is calculated as described above, but 
three additional parameters enter the calculation because 
the rapidity relaxation times and hence, the mean values 
of the fragmentation distributions as well as the widths 
of all three partial distributions and the particle numbers 
in the fragmentation sources differ from each other. If one 
takes the number of produced charged hadrons in the frag¬ 
mentation sources to be proportional to the corresponding 
numbers of participants as proposed in m, one needs only 
two additional parameters for asymmetric systems. 

Since particle identification is not available for pFh 
hadron rapidity distributions, the calculated rapidity den¬ 
sity distribution functions dN/dy have to be converted to 
pseudorapidity space, rj = —ln[tan(^/2)], before they can 
be compared to data so that the RDM-parameters can 
be deduced in minimizations. The Jacobian dy/dp = 
cosh(?7)[l + (m/pT)^+sinh^(77)]“^/^ that mediates the con¬ 
version depends not only on the masses of the particles, 
but also on the transverse momentum, and a careful con¬ 
sideration of both as detailed in na is required. The effect 
of the transformation is most pronounced at small trans¬ 
verse momenta, and in the mid rapidity range. 

In particular, it is not sufficient to consider only the 
mean transverse momentum, but rather precise results can 
be obtained by using the pion mass and calculating an 
effective transverse momentum {p^^) such that the exper¬ 
imentally determined Jacobian Jy=o of tho total charged- 
hadron distribution at ^ = 0 is exactly reproduced. This 

yields (p'"J^){b) = m^Jy=o/^l - 


with {pt) from the transverse momentum distributions 
for each particle species. For PbPb at 2.76 TeV we had 
obtained in 0-5% central collisions {p^^) — 0.323(pt)- 
The influence of the Jacobian on dN/dp at LHC energies 
remains, however, essentially confined to the midrapidty 
source, its effect on the fragmentation sources is marginal. 

Whereas the Jacobian produces a mid-rapidity mini¬ 
mum in charged-hadron pseudorapidity distributions for 
symmetric systems at LHC energies, its effect is not suf¬ 
ficient to generate the dip seen in the PbPb data at 2.76 
TeV, which is indeed more pronounced than the one in 
AuAu at RHIC energies - although the contribution due 
to the Jacobian decreases with increasing energy. It has 
been shown in |6] that an understanding of the mid-rapidity 
minimum requires both, the proper consideration of the 
Jacobian, plus the interplay of the three sources. In par¬ 
ticular, the fragmentation sources move further apart at 
higher energies, such that the gluon-gluon source deter¬ 
mines the mid-rapidity yield almost exclusively. 

In asymmetric systems like pFh at 5.02 TeV, the effect 
of the Jacobian is by itself not strong enough to produce 
a local minimum in the mid-rapidity source, see Fig. [T] for 
the minimum-bias pseudorapidity distribution of charged 
hadrons. Here we have used the same relation between 
{p^^) and {pt) as for PbPb. The Jacobian causes a de¬ 
formation of the mid rapidity source, and the local min¬ 
imum seen in data emerges when the gluon-gluon source 
is added incoherently to the two asymmetric fragmenta¬ 
tion sources to yield the measured pseudorapidity density 
distribution function. 

It is interesting to observe that the fragmentation sources 
become directly visible in measurements (SPS, RHIC) or 
predictions (LHC) of net-baryon or net-proton (proton mi¬ 
nus antiproton) rapidity distributions [20] , where the mid¬ 
rapidity gluon-gluon source cancels out. In a recent calcu¬ 
lation of the net-proton fragmentation sources that builds 
upon and extends the QCD-based model developed in [20] . 
results for 5.02 TeV pFh have been obtained [21]. Whereas 
the net-proton yields are probably too small to be mea- 


The effective transverse momenta are considerably smalleisurable, the charged-hadron yields from these sources are 
than the mean transverse momenta determined from the evidently substantial, and decisive for a proper under- 
Pt— distributions (Table [T]), and the corresponding effect standing of the pseudorapidity distributions of produced 
on the Jacobian is larger than what would be obtained charged hadrons. 
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3 Results compared with LHC data 

For minimum-bias (0-100%) pPb collisions at a c.m. en¬ 
ergy of y 5 AT AT = 5.02 TeV, the ALICE collaboration has 
already published data [3] , which have been analysed within 
the RDM in [ 6 ], see also the lower frame in Fig. [H Mean¬ 
while the ATLAS collaboration has also presented prelim¬ 
inary 0-90% minimum-bias data in [4], which we analyse 
accordingly in the upper frame of Fig. [H Here we have 
scaled the measured 0-90% ATLAS results by the geomet¬ 
rical cross-section ratio cr(0 — 90%)/cr(0 — 100%) ^ 0.898 
to make them directly comparable to the minimum-bias 
ALICE results. 



Fig. 2. (Color online) The predicted RDM pseudorapidity den¬ 
sity distribution functions for charged hadrons in minimum- 
bias pPb collisions at the future LHC c.m. energy of 8.16 TeV 
(upper solid curve) is shown together with the result at 5.02 
TeV (middle curve with ALICE data [3], see [6] and lower frame 
in Fig. [T]), and the minimum-bias distribution function in dAu 
at 0.2 TeV [22] with the RDM result from [17], lower curve. 


As is obvious from the figure, the ATLAS and ALICE 
results are compatible with each other within the system¬ 
atic error bars. Using the current version of the object- 
oriented data analysis framework ROOT [18], the RDM- 
parameters are determined through a minimization 
of the difference between the analytical solutions of Eq. [T] 
and measured dN/drj distributions. This yields, however, 
a larger mid-rapidity source for the preliminary ATLAS 
results (49 % of the total particle content), than for the 
ALICE data (35%). In view of the rather restricted range 
of the data in 77 —space, we expect that these discrepancies 
vanish once more precise final data in a larger 77 —range be¬ 
come available. 

With these forthcoming data, it will be very interest¬ 
ing to compare the predicted steeper slope on the proton¬ 
like side with the less steep one on the Pb-like side. A 
qualitatively similar behaviour had been observed in dAu 
collisions at RHIC energies of 200 GeV by PHOBOS [T]. 

At RHIC energy the mid-rapidity source in minimum- 
bias dAu carried only 23% of the total particle content 



Fig. 3. (Color online) The centrality-dependent RDM pseu¬ 
dorapidity distribution functions for charged hadrons in pPb 
collisions at LHC c.m. energy of 5.02 TeV are adjusted in 
the mid-rapidity region through minimizations to the pre¬ 
liminary ATLAS (left) [3] and ALICE (right) [5] data (sys¬ 
tematic error bars only for ALICE). The underlying distri¬ 
butions in the three-sources RDM are also shown, with the 
dashed curves arising from gluon-gluon collisions, the dash- 
dotted curves from valence quark-gluon events in the Pb-like 
region, and the dotted curves in the proton-like region (frag¬ 
mentation sources). 


m- This is in accordance with the expectation that the 
gluon-gluon moving particle production source that is al¬ 
most in statistical equilibrium, and centered at low rapidi¬ 
ties yeq in the fireball, becomes more important at LHC 
energies. 

Future LHC data on asymmetric systems will include 
charged-hadron production in pPb collisions at a c.m. en¬ 
ergy of ZiZ 2 l{AiA 2 ) • 2pp = 8.16 TeV for a 

proton momentum of Pp = 6.5 TeV/c. An RDM predic¬ 
tion for minimum-bias results at this energy is shown in 
Fig. [ 2 ] Here the particle content in the three sources, the 
corresponding mean values and widths of the distribution 
functions have been extrapolated based on the low-energy 
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minimum-bias values for dAu at 200 GeV min], aud the 
results for pFh at 5.02 TeV [6]. 

Siuce the euergy gap betweeu these two cases is very 
large, the extrapolatiou remaius somewhat uucertaiu. lu 
particular, the particle couteut aud the widths may turu 
out to differ slightly from what is showu iu Fig [3 where 
the iucrease iu the total yield - iutegrated over pseudo¬ 
rapidity - is rather small: The total charged-hadrou cou¬ 
teut iu miuimum-bias collisious is Nfot = 210 at 8.16 TeV, 
compared to 181 at 5.02 TeV. Note, however, that the dif- 
fereuce will be more prououuced iu ceutral collisious. 

Of particular iuterest is the iuvestigatiou of the ceutral- 
ity depeudeuce, which has beeu carried out experimeutally 
by both the ATLAS |1] aud ALICE [5] collaboratious. Sev¬ 
eral ceutrality bius are ideutical aud cau be compared di¬ 
rectly, with the followiug exceptious: For 0-5% the ATLAS 
values are deduced from the 0-1% aud 1-5% prelimiuary 
data m by weightiug with the geometrical cross sectious, 
similarly for 20-40%, which is deduced from the 20-30% 
aud 30-40% prelimiuary ATLAS data (here the geometri¬ 
cal cross sectious are ideutical for the two bius). No AT¬ 
LAS data are available for 90-100%. For ALICE, prelim¬ 
iuary data for various differeut ceutrality estimators are 
available. Here we have choseu the results from the es¬ 
timator that orders the eveuts based ou the uumber of 
clusters measured iu the secoud layer of the Silicou Pixel 
Detector (CLl, [5]). 

The results from the optimizatiou of the RDM- 
parameters for both ceutrality-depeudeut data sets are 
showu iu Fig [3 aud Table [T] (for ALICE). As au ad- 
ditioual coustraiut to the optimisatiou procedure we re¬ 
quire the RDM-parameters (meau values < yi ^2 > aud 
widths ri^ 2 ,gg) to bc mouotouic fuuctious of ceutrality, by 
choosiug appropriate iuitial couditious. The cousisteucy 
betweeu the two prelimiuary data sets, aud the corre- 
spoudiug RDM parameters is quite satisfactory. A small 
problem arises iu the 60-90% ceutrality biu, where the AT¬ 
LAS p-like maximum is above the Pb-like yield, but this 
deviatiou is withiu the systematic error bars. 

At all ceutralities, the mid-rapidity gluou-gluou source 
(dashed curves iu Fig[3|) has a larger particle couteut thau 
the sum of the two fragmeutatiou sources, with a rel¬ 
ative charged-hadrou couteut risiug mouotouically from 
58% iu the 0-5% ceutrality biu (referriug to the ALICE 
prelimiuary data) to 78% iu the 60-80% ceutrality biu. Iu 
coutrast, at a RHIC euergy of 0.2 TeV, the mid-rapidity 
source iu dAu had beeu fouud to have a larger particle 
couteut thau the sum of the fragmeutatiou sources ouly 
for 60-80% ceutrality [13 aud more peripheral collisious, 
where the valeuce-quark couteut is small aud heuce, the 
relative coutributiou due to the fragmeutatiou sources less 
importaut. 

For all three sources, the particle couteuts decrease 
mouotouically towards more peripheral collisious, whereas 
the fluctuatious, aud the associated values of the full widths 
at half maximum (8 Iu 2 )cr 2 ,^^ iucrease mouo¬ 

touically. The width Ti of the protou-like source decreases 
slightly for more peripheral collisious. These treuds are es- 



Fig. 4. (Color online)Produced charged particles in central 
dAu collisions at = 200 GeV (RHIC/ PHOBOS data 

[T]; see also [13) and 0-5 % pPb at y'SpN = 5.02 TeV (LHC/ 
ALICE data ID)- The upper solid line is a power-law fit to 
the data with the RHIC 0-20 % results scaled to 0-5% ac¬ 
cording to the cross section, and Ntot oc The number of 

charged particles produced in the fragmentation sources (trian¬ 
gles and dash-dotted curve) is Nqg oc In(spAr), and the number 
of charged hadrons produced in the mid rapidity gluon-gluon 
source (squares and dashed curve) is Ngg oc In^(spAr) (see text 
for details). The value of the pseudorapidity density for 0-5% 
central collisions at the equilibrium value of ij is also indicated, 
open squares and dotted line. 


seutially in line with what we had found earlier for dAu 
at RHIC energies m- 

For the centrality-dependent mean values of the frag¬ 
mentation sources, we have listed in Tabled] (^ 1 ^ 2 ) rather 
than the RDM-parameters since these can be cal¬ 

culated directly from Eq. Whereas the centrality de¬ 
pendence of the mean value of the proton-like distribution 
function is rather weak, a considerable drift in the mean 
value of the Pb-like distribution function is observed when 
the collisions become more central, and the interaction 
time tint increases. 

The mid-rapidity distribution is always centered at 
the impact-parameter dependent calculated values of i/eq, 
which are also listed in Table HI Obviously the centrality- 
dependent shapes of the total distribution functions are 
very sensitive to the precise values of the drift and also of 
yeq, which are properly described in the nonequilibrium- 
statistical RDM with a linearised drift term. 

Of particular interest is the evolution of the particle 
content in the three sources as function of the c.m. en¬ 
ergy. For the PbPb and AuAu systems, this has been done 
in [6] based on four RHIC central-collision data points 
between 19.6 and 200 GeV, and the LHC point at 2.76 
TeV. It has been found that the total charged-hadrou 
production is following a power law Ntot oc (satat/so)^’^^, 
whereas the particle content in the fragmentation sources 
is Nqg oc In (satat/^o). The particle content in the mid- 
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rapidity source is not too far from a power law in the 
intermediate energy range 0.1-2.76 TeV [6], but closer in¬ 
spection reveals that it obeys Ngg oc In^ (snn/sq) [H]. 

For a corresponding investigation in the asymmetric 
case, we presently have only two values, namely 200 GeV 
dAu and 5.02 TeV pPb. A double-logarithmic plot for the 
particle content in 0-5% central collisions is shown in FigS) 
The total particle content is shown as a power-law, Ntot = 

4.9- (5pAr/so)^‘^^ with sq = 10^ GeV^. The particle content 
in both fragmentation sources is Nqg = 26ln{spN/ sq) 
with So = 10"^ GeV^. The charged-hadron content in the 
gluon-gluon source is Ngg = 0.34 In^ {spN/ sq) with sq = 

1.9- 10^ GeV^. To assess how accurate these dependencies 
are, more data at different energies would be needed. 


4 Conclusions 

The analysis of charged-hadron pseudorapidity distribu¬ 
tions in 5.02 TeV pPb collisions in a nonequilibrium-statistical 
three-sources relativistic diffusion model reveals the charged- 
hadron content in the fragmentation- and mid-rapidity 
sources for particle production. We have determined the 
RDM-parameters in x^~ minimizations with respect to 
preliminary ALICE and ATLAS data. At all centralities, 
the mid-rapidity source has the largest particle content, 
but the fragmentation sources are necessary for a detailed 
understanding of the centrality-dependent shape of the to¬ 
tal distribution functions. It is the interplay of the three 
sub-distributions together with the effect of the Jacobian 
transformation that determines the pseudorapidity den¬ 
sity distribution of produced charged hadrons. For minimum- 
bias collisions, we have performed a prediction for 8.16 
TeV pFh that can be tested in the forthcoming LHC ex¬ 
periment. 

The shapes of the total distribution functions indicate 
that the system has not reached statistical equilibrium. In 
particular, the centres of the fragmentation distributions 
remain far from the equilibrium values yeq{b): It is not just 
the produced-particle yields that determine whether the 
emitting source is in statistical equilibrium, but rather the 
shapes of the distribution functions. A clarifying analogy 
may be found in the distribution function of the cosmic 
microwave background radiation, which - apart from the 
fluctuations - is indeed in perfect thermal equilibrium. 
This is, however, not the case for particle production in 
relativistic heavy-ion collisions when the full phase space 
is taken into consideration - although local equilibrium in 
the hydrodynamic sense appears to be achieved for events 
with pt < 6 — 8 GeV/c as may be inferred from the success 
of hydrodynamics to describe bulk properties. 
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